Armchair graphene nanoribbons: fP'T-symmetry breaking 
and exceptional points without dissipation 
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We consider a single-layer graphene nanoribbon with armchair edges and with a longitudinally constant 
external potential, pointing out that it can be described by means of an effective non-Hermitian Hamiltonian. 
We show that this system has some features typical of dissipative systems, namely the presence of exceptional 
points and of fPT-symmetry breaking, although it is not dissipative. 
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Since its isolation in 2004 fll], graphene has attracted a sig- 
nificant interest in the condensed matter community because 
of its unique electronic properties, the most astonishing one 
being the presence in its low-energy spectrum of two massless 
Dirac modes. These massless fermionic excitations allowed 
the experimental observation of exotic phenomena, such as 
the Klein paradox and the Zitterbewegung (see e.g. for a 
review), theoretically studied in the context of quantum elec- 
trodynamics long before. 

In this letter we analyze some properties of single-layer 
graphene nanoribbons with armchair edges in an external po- 
tential, pointing out a connection with the recently devel- 
oped theory of !P'T-symmetric non-Hermitian Hamiltonians 
in quantum mechanics. 

The systematic study of !P'T;symmetric Hamiltonians was 
initiated by the seminal paper [3], in which analytical and nu- 
merical hints were presented to explain the reality of the spec- 
tra of some non-Hermitian Hamiltonians. If the (P'T symme- 
try is realized in the spectrum, i.e. if every eigenstate of the 
Hamiltonian is also an eigenstate of the (P'T operator, it is 
simple to show that the energy spectrum has to be real. How- 
ever, since ^PT is an anti-linear operator, the symmetry can be 
spontaneously broken (for some enlightening examples and a 
review of the main results obtained in this field see e.g. 
and EPs (exceptional points) can appear Exceptional points, 
i.e. points for which two (or more) eigenvalues coalesce and 
the Hamiltonian is non-diagonalizable jstl, are a typical fea- 
ture of non-Hermitian Hamiltonians (see e.g. the reviews ^) 
with no Hermitian counterpart and they have been shown to 
produce experimentally observable effects iItI-B. 

To our knowledge, all the previously proposed physi- 
cal examples of systems governed by non-Hermitian fPCT- 
symmetric Hamiltonians involve dissipative systems, the main 
emphasis being on microwave cavities optical lattices 
ilCl lull and lasers il2ll . We show that, because of the 



spinorial nature of the wave function, some properties of 
graphene nanoiibbons can be described by means of an ef- 
fective non-Hermitian !f"T-symmetric Hamiltonian, although 
there is no dissipation. We give numerical evidence for the 
fPT-symmetry breaking and provide an order parameter Fi- 
nally we study the behavior of eigenmodes and eigenfunctions 
in the neighborhood of exceptional points. 



The model and the notations. We consider a nanoribbon 
section with open boundary conditions along the longitudi- 
nal X direction, armchair edges, and dimer lines across its 
width. The transverse distance (along the y direction) between 
the first line of lattice points not occupied by carbon atoms at 
the bottom edge and the analogous one at the upper edge of 
the ribbon {i.e. the "effective width" of the ribbon) is equal 
to L = {N + l)a/2, with a the graphene lattice constant. The 
wave function is equal to 

W) = Y.'^a{Ra)W-Ra) + iY,^B{RB)^{r-RB) , (1) 

where the (p(r)'s are the orthonormalized atomic orbitals 
of carbon, Ra and Rb are the positions of the atoms of the two 
sublattices A and B of graphene, and 
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(2) 



with K = 47t/(3fl) and j^, the 2-spinors corresponding 
to the two inequivalent Dirac points K and K' . The exter- 
nal potential V{x,y) is assumed to vary only in the transverse 
direction. Aside from its intrinsic interest, evaluation of the 
eigenfunctions and eigenvalues for such a potential represents 
the first step in conductance calculations for a generic poten- 
tial (varying in both spatial directions) performed by means 
of scattering matrix methods. Most existing calculations of 
transport in graphene based on such methods deal instead with 



a constant transverse potential ||13l [14|] 



We can expand in plane waves along the x direction 

and write = e^'^''^G;i'^ {y) and = e'^'-^(p ^'(y) . 

e.g. mm) 



The Dirac equation can be written in the form (see e.j 



r(/(y) + o,fc,-/o,8v)(p^(y)=0 



oA- + 'Ovav)9*^'(.y) = 



(3) 



where we introduced the shorthand f{y) = \V{y) — E]/\'f, 
with Vf the Fermi velocity (we use h~ \ in order to simplify 
the notations), and the armchair boundary conditions read 



(p^ (0) = 9"' (0) , cp^ (L) = e^iAL^^: 



(4) 
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The values for which this system admits non trivial solu- 
tions are the longitudinal momenta allowed in the nanoribbon. 
The problem of Eqs. (O-dlli can be rewritten in a more conve- 
nient form introducing, for y E [—L,L], the 2-spinor 



(5) 



and defining f{y) = f{L - \y\) and = fiy)ay - ikj,a,. 

We obtain 



idy^(j)=H'^{y)^{y), (p(L) 



2iKL 



9(-L) 



(6) 



which is formally equivalent to a Schrodinger equation for the 
non-Hermitian Hamiltonian VfH''\ if we interpret y/vf as the 
time. 

Symmetries. From Eq. (|6]l we deduce a simple result on 
the degeneration of the kx modes. If we denote by U{y) 
the time evolution operator associated to a given eigenvalue 
kx, i.e. (^{y) = U{y)(^{—L) with (^{y) in the correspond- 
ing eigenspace, the boundary condition can be written as 
t/(L)(p(— L) = e^'^^(\i{—L). By using the explicit form of //*^' 
it is easy to check that ' (L) = <3xU{L)<3x, from which it fol- 
lows that U{L)Ox^{-L) = e-^-'^^(5x^i{-L). If exp(4//:L) ^ 1, 
there cannot be eigenvectors of U{L) other than 9(— L) and 
Ox(p(— L), and hence just one independent eigenmode corre- 
sponds to each kx- From now on we consider lengths L for 
which this condition is verified, i.e. nanoribbons that are semi- 
conducting in the absence of an external potential also when 
edge relaxation is neglected. We denote by (p^.^ [y) the eigen- 
mode associated to kx. From the relations 

9A'J {y) - Ox {<^kA-y)) * (y) - (^z (9a, i-y)) * (7) 

it follows that if kx is in the spectrum then there are also k*, 
—k*x and —A:,. Thus the spectrum has a Z2 x Z2 symmetry. 

To reveal the ^'T symmetry of this problem it is convenient 
to take the square of Eq. (|6]l and project 9 on the eigenstates 
of <3y. If we denote these projections by they satisfy the 
equations 



Pl - p T ^± (y) = -kl^± (y) : 



(8) 



which are clearly invariant under the T'T transformation, be- 
ing the action of the operators and 'T defined by py — > —py, 
y —J- —y and py ^Py, y ^ y, i ^ ^i, respectively. If the 
(P'T symmetry is unbroken then kx has to be real or imaginary; 
complex conjugate (intended here as a number with nonzero 
real and imaginary part) pairs appear in the spectrum only if 
this symmetry is broken. 

We explicitly notice that if the Schrodinger equation is 
used instead of the Dirac one, the equation corresponding to 
Eqs. (|3]l is similar to Eq. (|8]l but with an Hermitian left hand 
side, so that all the kx values have to be real or imaginary. 

In the presence of spontaneously broken symmetries it is 
customary to look for an order parameter, i.e. an observable 
that vanishes when the symmetry is realized in the spectrum. 
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Figure 1. Coalescing of two eigenvalues along the real line for 
e = 1/5; the arrows denote the direction of increasing a values 
(Aa = a — 0.060953). In the inset (a) the imaginary parts of the 
coalescing eigenvalues are shown together with a fit of the form 
a = ci + C2{lva{kxa))^. The inset (b) shows the behavior of the 
square of the transverse momentum together with a linear fit. 



a non zero value signaling the symmetry breaking. We point 
out that the mean value of the transverse momentum 



Py^ ^>{yy{-id,)^>{y)Ay 



'P(y)T'P(y)d3; (9) 



satisfies this requirement. This can be easily proved exploit- 
ing the symmetries of our problem; from Eqs. (|7) and be- 
cause [^(y)^ is !P-even and (p(— 3')"^Oz9(3') is constant, it fol- 
lows that the numerator of Eq. (|9]l vanishes if k^ is real. We 
checked numerically that the transverse momentum py is dif- 
ferent from zero when kx is complex {py appears to vanish 
also for complex kx values only when the potential satisfies 
V{y) ~ V{L~y)). Thus, in this system the realization of 
the symmetry in the spectrum is related to the value of 
the transverse momentum, which is an order parameter for 
the fPCT symmetry breaking. This is completely analogous to 
what happens in QCD with compact dimensions, when charge 
conjugation can get spontaneously broken and the related or- 
der parameter is the baryon current in the compactified direc- 
tion ifieli . 

Finally, it is worth mentioning that the transformation that 
reverses the sign of f{y) {i.e. <5^) is unitary and indepen- 
dent of Py, and hence |9(3')P is invariant for simultaneous 
flipping of the signs of the potential and of the energy. This 
is just a manifestation of the chiral symmetry of the system. 
As long as /(y) does not have a definite sign, heuristic argu- 
ments based on that symmetry indicate the presence of eigen- 
functions localized around the minima as well as eigenfunc- 
tions localized around the maxima of the potential; the former 
(latter) ones are expected to describe particles with positive 
(negative, i.e. of opposite sign with respect to ^v) group ve- 
locities. These observations provide a simple argument for 
the existence of some singular behavior: by varying the en- 
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ergy, eigenmodes with opposite group velocities coalesce in 
a non-analytic way, because they cannot combine in a two- 
dimensional eigenspace. 

'2^ -symmetry breaking and exceptional points If f(y) is 
constant, Eq. (|6]l can be analytically solved: all the values 
are real or imaginary and all the eigenstates of the Hamilto- 
nian (their projections on the eigenstates of o, , to be precise) 
are also eigenstates of the operator This is no longer true 
when an external potential with non-trivial y dependence is 
present and we now provide numerical evidence that in this 
system !P'T-symmetry can get spontaneously broken. As a 
simple example we use the Lorentzian shaped potential 

125 

= 250 + a2((v/«)- 150)2' 

with a,e e M. The general qualitative features are however 
independent of the particular potential adopted. We choose 
L = 500fl as the effective width of the nanoribbon. 

It is simple to (numerically) check that by varying the pa- 
rameters (a and e in our example) two different phenomena 
can occur: 

A) the number of the real eigenvalues varies but the number 
of the complex ones is preserved; 

B) the number of the complex eigenvalues changes. 

In case A) a couple of real eigenvalues turns into a couple of 
imaginary ones (or vice versa); the Z? x Z2 symmetry then im- 
plies the existence of some values of the parameters for which 
^jt = is a doubly degenerate eigenvalue. However, we are 
assuming exp(4;/rL) 7^ 1, and hence only one independent 
eigenfunction is associated to each eigenvalue; as a conse- 
quence, this point of the parameter space is an exceptional 
point. Case B) is completely analogous: two real or imagi- 
nary eigenvalues coalesce and become complex (see Fig. [T]i. 
However, while in case A) the ^^'T symmetry is unbroken (or 
broken, if the number of complex eigenvalues is different from 
0) irrespective of the EP, in case B) the EP is associated to fPCT 
symmetry breaking. Clearly this kind of EP appears in a spec- 
ular way both in the upper (right) and lower (left) half-plane, 
because of the Z2 x Z2 symmetry. 

Behavior near EPs. If we are not at an exceptional point, 
the variations of the eigenvalues and eigenvectors are linear in 
the variation 8/ of /; in particular it can be shown that 

J_^dy^li-y)-^^^-My) = - J _^dy^i{-y)^kAy) ■ (H) 

The group velocity vf in the x direction is obtained in the spe- 
cial case of an energy shift and is given by 

8£ J-Ldy^lX-y)<?kAy) 

t>k^ /_Ld3'<.(-}')av(pA,(3') 

The velocity defined as above has a clear physical interpreta- 
tion for real modes, for which it is real; however also complex 



group velocities can provide interesting information in many 
physical systems (see e.g. iHtIi ). 

At an EP, /f^dytp^ (— y)^^^. (y) vanishes for some kx and in 
the neighborhood of an EP corresponding to the eigenvalue 
we find 

[kx-k'^x'^f ^ 1 S\Ay^l{-y)Ox^kAy) 

' ^ 'Fi:.I-LdyV,S-ynkAy)' 

(13) 

This is nothing but the well known square root behavior, in the 
neighborhood of an EP IH^], of two coalescing eigenvalues 
as a function of the external parameters (see the inset (a) of 
Fig. [T]i. Near an EP we can then factorize the Hilbert space 
as the product of the 2d space of the collapsing eigenfunc- 
tions, for which the rapidly change with a small shift in E 
(fast modes), and the span of the other modes (slow modes), 
which can be assumed as fixed in the neighbourhood of the 
EP. Here we are assuming that eigenfunctions associated to 
different exceptional points do not mix; in the considered nu- 
merical examples we checked that this assumption is indeed 
true. Thus: 

where klc^\ k^x'' are the two coalescing eigenvalues, 7? is a 2 x 2 
matrix and 9± are the two initial states. Notice that if the 
difference between one or more eigenvalues and the coalesc- 
ing ones is 8^v ^ /J, they can mix together It turns out that 
if they are quasi-degenerate with the coalescing eigenmodes 
(bkx <C /j) the mixing between them is just a rotation, that is 
irrelevant for the features that we are going to describe. As 
long as we are not at the EP, we can choose the normaliza- 
tion /_^^d(p(p^, (— 3')(p*:,.(3') = ?>k^k[^ so that RR^ ~ I. For the 
sake of simplicity, in the following we restrict to the case of 
8f = —dE /vp and EPs on the real axis. The qualitative results 
obtained are nevertheless of general validity. Before the EP 
is reached, the group velocities of the coalescing real eigen- 
modes are opposite in sign, as shown by Eq. (fTSl l and indicated 
by the subscript ± in Eq. ( fT4b . and from Eqs. d?) it follows 
that R* = 0-R0-. After the EP is crossed, the eigenvalues as 
well as 5E /5kx become complex conjugate and from Eqs. (|7]i 
it follows that R* = ztOxRo-. Relaxing the normalization con- 
dition, we can parametrize R as follows: 

R^l + e'%.. (15) 

The domain of definition of the parameter 9 is ]°oi, Oi] U [0, 7i] U 
[n + Oi,n + °o! [. The kx values are real if Im 9 is different from 
and complex if 9 is real; the larger Im(9) the further apart 
the modes kx ' are. When 9 = the eigenfunctions (p (i 21 are 
linearly dependent, so this value corresponds to an exceptional 
point. We rewrite Eq. (fTjt in terms of the parameter 9 in the 
simplest case in which the term A = /^^ dy(^+ (— y)-^av(p- (y). 
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which appears in /v, is neghgible 



1,2) 



'V_')cos9=F'(i'+' — v_')sin9 



(v;'-v:') 



-h2 



COS 9, 



(16) 

where and £0 are constants and v± are the group velocities 
associated to the eigenfunctions (p±. The condition A w is 
found for example when one eigenfunction is localized around 
the minima and the other around the maxima of the potential. 
The approximate Eq. ( fT6] l is accurate only in a neighborhood 
of the EP with 9 ~ but, if /j is much less than the energy 
scale in which the Hilbert space factorization remains valid, 
then the value 9 = 7t of Eq. (fTsT i corresponds actually to an- 
other EP, and the previous approximation is good in the whole 
interval < 9 < 7t. This happens when two real eigenval- 
ues collide, become complex and then come back on the real 
axis, the two EPs being sufficiently close to each other In this 
case Eq. ( fTSI l captures the whole out-of-axes "motion" of the 
eigenvalues. In Fig. [2] the results predicted by Eqs. (fTST i and 
(fT6] l are checked against numerical data and the agreement ap- 
pears to be excellent. After crossing both exceptional points, 
the eigenfunctions almost return to the starting ones; observe, 
moreover, that the energy scale of the phenomenon shown in 
Fig. |2]is of order sa IQ^^vp/a, to be compared with an an- 
alytical background of order w Q.lvp/a (see the caption of 
Fig|2]i. These aspects make the numerical observation of the 
phenomenon extremely difficult, so that the effect of two very 
close EPs may be incorrectly interpreted as a mode-crossing. 

Conclusions. We have pointed out a connection between 
properties of graphene and the theory of non-Hermitian 
Hamiltonians, by showing that armchair graphene nanorib- 
bons provide the first example of a nondissipative system de- 
scribed by a !PT-symmetric non-Hermitian Hamiltonian. We 
have also established that the transverse momentum py is an 
order parameter for the symmetry breaking. We have 
numerically verified the presence of exceptional points and 
shown that, in their neighborhood, the behavior of and of 
the eigenfunctions of the Dirac Hamiltonian, in the presence 
of a longitudinally invariant external potential, is theoretically 
well understood. 

An aspect that certainly deserves further study is the effect 
of exceptional points on the transport properties of graphene 
nanoribbons, in the presence of a potential that varies also in 
the longitudinal direction. Moreover it would be interesting 
to study more in depth the properties of complex eigenmodes 
and the effects of the non-vanishing transverse momentum. 

We thank M. D'Elia for useful comments. D. L., M. M., and 
P. M. gratefully acknowledge support from the EU FP7 1ST 
Project GRAND (contract number 215752) via the lUNET 
consortium. 
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Figure 2. The momenta (top) Aytj''^' ^ 4''^' - 0.199894795/a and 
the energy {bottom) AE ~ £ + 0.20065 1930 /a of two coalescing 
modes as functions of cos 9, estimated from (deti? — 2)/\/l — delR 
(see Eqs. ( il4b and dlSb). for the potential JlOt with a = Jl/2 and 
e = —aE/vf. The lines in the upper graph follow the prediction 
l il6t : the asymptotic velocities v± are obtained by considering AE > 
2 ■ 10^^ Vf/a and jj obtained from the linear fit of the energy versus 
cos 9 shown in the bottom graph. 
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